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ABSTRACT:  An  investigation  of  the  basic  relationships  which  describe 
the  flow  of  non-Nevtonian  fluids  under  conditions  of  impact  loading, 
possible  plastic  flow,  and  short  periods  for  the  various  conditions. 


A  generalization  of  the  Navier-Stokes  equation  is  proposed 
for  non-Newtonian  fluids  in  which  the  stress  is  proportional  to  a  power 
n  of  the  rate  of  strain.  Results  show  that  the  ratio  of  the  transverse 
to  the  longitudinal  3traln  is  not  a  constant  but  a  function  of  time  or 
of  strain  and  depends  rtrongly  upon  the  visco-elastic  response  of  the 
medium  as  well  as  on  the  loading  conditions.  The  pronounced  effect  or 
volume-viscosity  on  ti  <•  shape  of  this  function  is  illustrated. 
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Pseudo  Plastic  and  Dilatant  Fluids 

flie  s Jjiple  s hear  flow  "of  "pseud o-nlastic  and  dilatant  fluids 

is  expressed  with  considerable  exactness  by  Waele-Ostwald's  formula 

as  follows:  .  n 

du  Tn 

^y"^pSU  (i) 

where  yUpSU  is  the  quantity  similar  to  the  viscosity  of  Newtonian 
fluid,  r  is  the  shearing  stress,  du/dy  is  the  velocity  gradient, 
n  is  the  rheological  constant,  7 he  flow  of  pseudp-plastic  fluids  is 
expressed  by  n  >  1,  and  that  of  dilatant  fluids  is  given  by  n  <1 
■  in  the  above  equation,  respectively. 

The  relation  between  the  stress  components  and  rate-of- 
s  train  components  for  non-Newtonian  fluids  can  be  expressed  as: 

Pik  *  ^ik  -  P^ik.  '  •;  (2)  . 

There  is  a  component  of  the  stress  tensor  ( r >  as.  indicated 
in  the  following  array: 


(Pik> 


and  eik  is  a 
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These  tensors  are  symmetric  tensors,  that  is: 

Pik’Pki*  eik-,ki  (1'lk> 

p,.  is  a  normal  stress  if  i  ■  k,*a  tangential  or  a  shear  stress  if 
*K  i  /  k,  p  is  the  pressure,  which  is  the  mean  of  the  normal  stresses 
over  three  planes  mutually  at  right  angles.  6^  is  the  Kronecker 
delta  (6^-  1  if  i*k,  6^k*  0  if  i-k ;  i,k»  x,y,  z).  u,  v,  w  are 
parallel  velocity  components  to  the  coordinate  axes  of  rectangular 
Cartesian  coordinates  x,y,  z,  respectively,  e  ,  e  .  etc.  ere 
Au  e  .  av  +  au  **  yy 

exx"  W  »  xy  W  3y“  etc.  (A  denotes  partial 

derivative) 

^is  an  arbitrary  scalar  function  of  the  three  invariants  of  the 
rate-of -strain,  and  may  be  expressed  as 


e  4  e  ♦  e. 

tY  '  W  .  1 


i  2  2  2 

e  •  e  e  +e  e  ♦  e  e  -*(e  +e  ♦  e_„ 
yy  zz  zz  xx  xx  yy  4V  yz  zv  xy 
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)  '*) 


|  e ..  ]  •  e  e  e  ,  e  e  e 
1  ik1  xx  yy  zz  ♦  X  yz  zx  sy 

t(  «  e2 

4  xx  yz 


♦  e  e 
yy  zx 


♦  e  e  c  )  (7) 

zz  xy  '  ' 


If  we  treat  the  fluid  as  being  incompressible,  the  first  invariant 
0  vanishes,  or  ^ 

9  0  A  a  A  A  v  M  l  <  I 


e  ♦  e  +  e 
xx  yy  zz 


Consequently,  T^can  be  expressed  as  an  arbitrary  scalar  function  cf 
the  two  invariants  e,  |e*.  |.  in  rectilinear  flow  and  in  two-dimen¬ 
sional  flow,  the  third  iiWariant  j e^.  |  is  identically  zero.  Hence 
^  can  be  expressed  as  a  scalar  function  of  the  second  invariant  e 
only  in  three-dimensional  flow.  Taking  account  of  eq.  (£),  the 
second  invariant  e  becomes 

e  ■  -  i[  2(e  ^  +  e  ^  +  e  +  e  ^  +  e  ^  ]  (?) 

4  xx  yy  zz'  yz  zx  xy  J 

If  we  generalize  ’Vaele-Ostwald’s  formula  by  assuming  that  the  stress 
components  are,  in  general,  expressible  by  the  l/n  th  rower  functions 
of  the  rate-of-s train  components  for  non-Newtonian  fluids,  the  rela¬ 
tions  between  the  stress  components  and  the  rate -of -strain  components 
may  be  expressed  as 

Pii  ■  -P  +  eii  (i-x,  y,  z) 

'  U  P  (10) 

^  0,  (i,k,  ■  x,  y,  z;  i  {  k) 


Pik  - 


2  2  2,  2 
K'svv  +  e,„t  +  e_„)  4  e  4  e 
XX  .  yv  ZZ  yz 


2  +  e  2  H2" 

fV  YV  J 


If  we  include  the  kinetic  reactions  and  the  external  forces,  and 
consider  the  equilibrium  of  an  infinitesimal  rectangular  parallele¬ 
piped  with  its  centroid  at  a  point  Pf'x,  y,  z,)  and  its  edges  parallel 
to  the  ’ordinate  axes, then  by  resolving  the  resultant  forces  along 
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the  axes,  we  obtain 

p  £“.  ■  (x  ♦  4p«c  ♦  *?»♦  *p*» 
r  Ot  Z5 —  Zy  IT" 


where  X,  Y,  Z  are  the  components  of  the  extefcnal  force  actin'*  on  the 
fluid  per  unit  mass,  Q  is  the  density,  and  Du/Dt,  Dv/Dt,  and  Dw/Dt 
are  the  components  of  the  acceleration  of  a  particle  of  the  fluid. 
Substituting  eq.  (10)  in  eq.  (12),  then  making  use  cf  e"  (c)  we  hrve 
the  fundamental  equations  of  motion  for  non-Newtonian  fluid  in  the 
form 


?  BT  *  P*  ■  H  [  °^u  4  2  e„j7  *  «Jy  *  •„  jr|] 

m  r .  .  /  r  _  2  iO  -  AQ  fiOt 

?rc*Pr-$*  #7V4 


•ftt 


where 


l  *  ezx  sl  +  ezy  27  *  ?eK  ET,[ 

2  •  2 
A  .  A 


u  .  w  .  w 

is2  4  ~s  *  ~s 


The  equation  of  continuity  is  evv  +  +  e„_  "  ®  (S’) 

xa  yy  zz 

This  is  treating  the  fluid  as  being  incompressible.  Parameters 
should  be  theorized  for  variation  in  density  for  fluid  under  pressure 
and  eqs.  (12)  and  (13)  should  be  solved  with  the  appropriate  terms 
inserted.  This  will  be  investigated  in  another  papei. 

If  n  »  1,  eq.  (13)  becomes  the  Navier-Stokes  equation  for  Newtonian 
fluid.  The  rheological  constant  n  has  the  different  values  for 
varied  fluids  and  n,  in  general,  seldom  has  an  integral  value.  The 


approximate  solution  by  the  energy  method  may  be  applied  for  non- 
Newtonian  flow.  The  dissipation  of  energy  in  non-Newtonian  flow  may 
be  calculated  by 


latea  oy 

1  "  J  J  J  ^exxpxx  +  ^TpTJ  *  ezzpzz  +  eyzpyz  +  ezxpzx+exypxy^ 

[dxdydz] 

Substituting  eq.  (10)  into  this  equation  and  then  using  eq  ('”)  we 
I"A,p^n)  jf  [2(exx  *  ew  +ezz)4eyi4e2x*e^*^ydl  Ul,) 


e  d  +#d  +ep  +ep  +ep+e 
xxpxx  yyyy  zz'zz  yzKyz  zxhzx  xj 


so  that  the  dissipation  function  0,  which  is  the  rate  cf  dissipation 
of  energy  per  unit  time  per  unit  volume  is  obtained  as 


v  1/n  r<9,  2  2  2x  2  2  2^+l>2n 

0  ■  li  '  [2(e  +e  +e)  +  e  +  e  +  e  T 

r  'psu  1  xx  yy  zz'  yz  zx  xyJ 


Now  if  the  velocity  distribution  of  flow  gets  the  minimum  dissipation 


HOWARD 


of  energy,  that  ie  61  ■  •///<  dxdyd2  "0  (16) 

It  may  be  proved  that  the  velocity  distribution  also  satisfies  the 
equation  of  motion  in  which  inertia  terms  are  neglected  and  external 
forces  derived  from  a  potential  or  zero,  that  is 

9  +  2exx  H  *  exy  A?  +  exz  AI  "  ET 

9  *  V  E  t2°yy%  *  'y*  ZT  ‘  §  (17> 

2  AQ 

Ortw+e  A©  0  AO  AH 

V  «  17  *  Sy  Jy  *  K  ’  E 

Where  H  is  a  function  of  x,  y,  z.  From  eq.  (15)  we  have  • 

(0  .  !Lti  u  Vn  g(j*“  M»  .2  “I  ♦2$*  *2!! .  ({S  ♦  |»}(|SS,  fiSlj 

n  'pr:  1  Ax  Ax  Ay  Ay  Az  Az  Ay  AzAy  £z 
+  (6u  +  AW)  /J6u  +  ASW)  4  4  JU)  (J6l  4  J6U)  ] 

'I!  ff  lET  ET'  'Ax  Ay'  'ST  AT'  J 

Therefore  the  variation  61  of  I  becomes 

at  .  n+1  «  1/n  f  f  fan  Au  A8u  Av  A6v  «  Aw  A6w  .,Aw.Ay, 

41  —  /"pm  J  J  J8[2  E  ET  4 .  2 -5?  Ty  *  2  15  iT  *(arE) 


(A6w  A6V)  (ou  aiiv  (cou  *  oow,  ,nv  aux  /flovaou) ) 

V  2x;  ^Tz  Ax'  'Ax  Ay;  vAT7r'J  y 

Considering  that  Gu,  6v,  6v  are  zero  at  the  body  surface  and  infinity 

and  then  using  the  integration  by  parts,  the  differential  terms  of 
6u,  6v,  and  6w  with  respect  to  x  become 


r6u  .  AWv  ,A8u  ♦  A6v 


Av  .  Au\  /A6v  A6ux 


j  j  j  ox  ox  ox  oy  ox  -  oz  ox  jj” 

a.  W-. 


+(z? 4  £r)6w3d3{  “i0^  *"  +'n(Ar+  :^6v'  (fe)(AivH)6M}^i;dydz . 
"J//  ^25uK  ‘  6v!H  *  $  *  H11  dxdl,dJ ; 

-  J  J7 8 126“  §  * {  n  *  ® .*»  • *  &»***..  •  •  ■  ■ 

Since  the  differential  terms  of  5u,  6v,  and  6w  wtih  resrect  to  y 
'  and  z  are  transformed  in  the  same  way,  eq. .IF  may  be  written  as  ‘ 

51  .-UkL)^'1/11  f  f  f[6u(0^2u  +  2e  ~  +  e  ~  +  e_  £)  ♦  6v(  O^v 
n  psu  JJJ1  v  xxAx  xy  Ay  *z  Az  v 


•  AO  .  AO  AO 

+  e  -r—  +2e  t'-  +  e  — - 

xy  Ax  yy  Ay  yz  Az 


AO  AQ\  e  2 

S*  ev,  *  ^(Sy-w  ♦  e 


AO  AO  0  /!0» 

t-  +  e  T-  +  2e  —  J 

zx  Ax  zy  Ay  zz  rz 
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Utilizing  eq.  (6*)  and  eq.  (16)  the  above  equation  becomes 
f  [6u(®  *  ?exxTx  *  8xy  I?  *  •»  if1  *  6V  (9  ff  ♦  2«w 

17  *  V  K>  *  6WM  ’  e2x  E  *  esy  17  *  2«„  ZT,]  4xd^J  "  0  M 

Since  the  variations  6u,  6v,  and  6v  must  be  satisfied  by  eq.  (6)  so 
that  Abu  A6v  Afiw  n 

ST  !y“  Is"  ’  0 

Multiplying  this  equation  by  a(x,  y,  z)  and  then  adding  its  resulti 
to  eq.  (19)  produces  the  following  equation!  (20) 

JJ J  IM«v!u  *  2exx  n  *  e«y  *  exx  I>  *  6W’  *  VS* 

Jejry  r7*VH''*MCy?w,a2x  sK,  h>‘  a(ir*i7!4rlta,lydl 

equals  zero.  Since  6u,  6v,  and  6w  are  zero  at  the  body  surface  and 
infinity,  using  integration  by  parts  gives 

I/J  a  jjbudxdydz  ■JJ[|a6ul  -J^budx]  dydz  ■  •i/I  ^budxdydz 

Therefore  eq.  (20)  becomes 

f  f  f((W  t^u  +2e  -p+e  ~  +e  -  t-)6u  +  (W  erv  +e  £  *2e  £ 

J  J  J  xx  Ax  xy  Ay  xz  Az  Ax7  v  yx  Ax  yy  Ay 


yx  Ax  yy  Ay 


+e  —  -  t^6v  +(0  rrw  *e  +e  '  42  +2e  42  -  febwjdxdydz  »  0 

yz  Az:  Ay  v  .  zx  Ax  zy  Ay  zz  Az  Az 

As  6u,  5v,  and  bw-may  take  arbitrary  values,  in  the  domain  of  flow, 
the  coefficients' cf  6uj  '6v,  and  bw  must  be  zero  in  order  to  always 
•  satisfy  the.  above  equation.  Then,  the  same' equation  as  eq.  17  is 
•obtained  as  fellows; 

*  o  2  '  o  AO  .  «  AO  .  •  AO  Aa 
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Thus,  the  velocity  distribution  that  satisfies  e  ].  (16)  is  found  to 
be  the  solution  of  eq.  (17),  if  the  inertia  terms  are  neglected  and 
external  forces  derived  frem  a  potential  or  zero.  Therefore,  if  the 
flow  satisfies  the  conditions  mentioned  above,  the  approximate  solu¬ 
tion  Of  flow  may  be  obtained  by  using  this  minimum  dissipation  of 
energy  method;  the  velocity  distribution  ccr  taining  some  undetermined 
coefficients  and  satisfying  the  boundary  conditions  is  api rot riately 
assumed,  and  then  the  undetermined  coefficients  are  det’rmined  so 
that  the  dissipation  energy  becomes  minimum. 
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VISCO-ELASTICITY 

Some  of  tho  characteristics  of  visco-elastic  substances  are 
desirable  in  the  operation  of  hydraulic  systems  of  aircraft,  weapon 
systems  and  transportation  equipment.  In  effect,  there  is  a  need  to 
absorb  energy  without  a  shock  result.  In  viscc-elastic  substances, 
the  ratio  of  lateral  to  longitudinal  strain  in  uni-axial  stressing  is 
not  a  constant  but  a  variable  parameter  decendinc  on  tine  and  the 
conditions  cf  the  test.  A  theoretical  analysis  cf  this  rarnmeuer  un¬ 
der  various  assumptions  concerning  the  visco-elastic  response  of  the 
material  and  the  testing  conditions  as  developed  by  Freuder.thal  and 
others  is  given.  The  Maxwell  body  mere  closelv  ar'-reximstes  the 
desired  reaction.  Promise  of  further  arnroxima.ien  is  shown  in  the 
characteristics  of  some  long  chain  polymeric  fluids  which  show  a  hi ih 
bulk  modulus  under  high  pressures,  i'he  final  solution  to  the  problem 
will  be  accomplished  by  development  of  the  desirable  characteristics 
from  a  molecular  structure  approach. 

The  stress  and  strain  components  o^  and  may  be  written 

as  S..,  e..  and  a.  ,  £*  ,  respectively,  so  that 

-  *«•■«’  *"  tl>  "•*' 
where  6..  is  the  sane  aronecker  delta.  The  general  linear  visco-elas¬ 
tic  medium  is  defined  by  the  linear  operator  equations 

F  •  20®  and  *0^  •*«■(£,*- 3ar)  (!•*) 

where  T  denotes  the  tenperature  difference  with  respect  to  a  constant 
reference  temperature  and  c  the  coefficient  of  thermal  expansion. 

The  constants  of  the  operators 

. . . 


a*  ♦  ai 


1  K* 


\7.ln 

*'rW 


(1.?) 


are  combinations  of  relaxation  times,  retardation  times  and  shear 
moduli. 

Combining  eqs.  (1.2)  under  the  assumption  cf  isothermal  condit¬ 
ions,  the  visco-elastic  stress-strain  relations  may  be  obtained  in 
the  altemativeforms: 

" (3K  F*20  rW3)£kit  8ij  •  20  r  ^t) 


and  . 


“““  '  c  ./IP  1  F  x  _  .  .  .  1  F  „ 

^ij  T  “  &  2  '  akk  6ij  +  2&  $  °ij 

Poisson1  s  ratio  for  the  elastic  medium  related  to  K  and  G  by  the 

relation  ^  .  {%  -  2G)/  2(3K  +  G)  (1.5) 

can  be  transformed  and  becomes,  in  linear  visco-elastic  theory,  the 


(U) 
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OPERATOR 


■  3KP  '5*  -  2G?P» 


(1.6) 


^  2(3K?  *  ♦  art  •)  v  ' 

Dividing  the  numerator  and  denominator  by  1?  KO,  ei.  1.6  takes  the 


1  P  _  1  IP 
) _ ^  ?  a i? 

I  ♦  i  f 


(1.6a) 


9);  ? 


Ti  T 


which  can  be  obtained  Ly  evaluating  the  second  eq.  (l..'i)  fcr  the  con¬ 
dition.  Thus  ,  .  « 


.  .  ( I-H.  IJL  ,  , 

'h  9K  ;■  3d  ;  1 


9K  ?  60  Q  11 


(1.7) 


the  time  dependent  ratio  >)(t)  is  obtained  as  the  ratio  too  ■ 

-  22/  t  u  *  V  ( CTll) 

Laplace  transforms  are  used  to  evaluate  this  ratio  by  use  of  the 
stress-strain  relaicns  (1.2)  and  (l.'i)  and  the  operator  equations  (1.3) 


T*  (p)  -  £.  a,  pk 
k-1  K 

k'r  k 

9*  (p)  m  pK 

k-1  K 


’’Mp)  •  £  b.  p^ 


*  'p) 


(1.P) 


The  ratios 


l  P  (r)  ,  *'  (?) 

—  (p)  -  *  /  x  and  (p)  - -  (1.9) 

B  "(r)  r  t*  (-p) 

are  quotients  of  two  polynomials,  the  inverse  transferns  of  which 
can.  usually  be  expr  ssed  in  the  form  of  series  of  negative  exponentials 


The  transformed  ear.  (1.7)  can  be  directly  written  in  the  fora 


£u  -  (- — -- — 
9K  K'  (p) 


%  T  (P) 


% 


33  a  (r) 


-  )  O. 


6G  K  (p). 


(l.](D) 
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Poisson's  ratio  >J(t)  is  obtained  as  the  inverse  transform 

_ 1  1 


ir\[ 


%  W  (p)  60  i  (p) 


j  v 


(1.11) 


L'A([ 


]  * 


11 


%  (p)  30  1  (p) 

When  £  jj  ■  £n  (t)  and  (p)  is  given,  with  -  0, 

is  evaluated  with  the  aic  cf  the  first  eq.  (1.10)  and  introduced 
into  the  second.  Hence  the  ratio 


(t) 


9  K  K»  (p)  6  G  1  (p) 


L  9  K  1*  (p)  3  0  1  (p)J 


e  (t)  d.i2) 

H 


Eq.  (1.11)  and  (1.12)  will  now  be  evaluated  for  the  Harwell  type  of 
idealized  linear  visco-elastic  behavior  in  volume-constant  distortion. 
If  we  assume  elastic  voluajetric  deformation  1*  reduces  to  unity.  The 
testing  conditions  investigated  are  :  constant  stress  <7^  ■  const., 

constant  stress-rate  ■  ct  where  c  is  an  rbitrary  constant, 

constant  strain  6^  •  const,  and  constant  strain-rate  »  ct* 

For  the  actual  evaluation  the  arbitrary  relation  3  K  »  p  0  has  been 
introduced,  for  which  the  frequently  usel  elastic  value  >)  ■  l/3 
is  obtained. 


Evaluation  of  "Fcisscn's  Ratio  forSlrle  Elastically 
Compressible  v i sco-Elas  ETc d  1 a^? fax  well  ^cdy 

The  mechanical  behavior  cf  a  Maxwell  body  is  determined 

by  *  -  \  +  P  5  *  P 

P'  -  1  '  Q'-  1 

where  t  ■  **[/  G 


Constant  stress  is  specifically 

an  *  c  H  (t)j  ^ 

From  eq.  (1.11) 


i 


* 


I 
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t,,  20 

.  1  T*1-^ 

r 

?  * 1  *  rr 

i’he  diagram  of  ^(t)  for  G/K  ■  3/3  is  shovm  in  figure  1,  curve 
*c.  At  t  ■  0,  the  value  of  ^  ■  1/3  and  as  t  0$  ^  ■  1/2. 

These  limits  may  be  explained  by  the  fact  that  the  instantaneous 
response  cf  a  Maxwell  body  is  purely  elastic  snd  consequently  is 
controlled  by  the  ratio  G/K  only.  Inder  long  time  loading  a  Max¬ 
well  body  behaves  as  a  viscous  fluid  which  results  0-*l/2,  with 
the  elastic  compressibility  becoming  unncticeable  under  the  uni¬ 
axial  testing  conditions. 

The  constant  stress  rate  is  specified  by 

°n  ■  ct  H(t)  and  ^  ■  c  -2 


Frcm  eq.  (1.11) 


L-1[( 


1  +  T 


6  G  t 


:>-r> 


l*1  [(Hi* . -L> -^] 

3  G  t  p  9  K  p* 


t  .  2  C- 
7T  1  ‘  IT 


2  +  1  *  JT  ) 


The  corresponding  diagram  is  shown  in  figure  1,  curve  ■  ct 
Constant  strain  is  specified  by  _ 

■  c  H(t)  and  "  C  p 

^rom  eq.  (1.12) 

.  ,,  3  K  (  i  *  p)  -  2  Gp 

-0  (t)  -  L  1  [ - i -  ] 

2[  3  «  i  *  p)  ♦  G  p  ]  p 

-  ?  ♦  ("y  rir~  - 1 )  «f  (-  jj-tg  "? }  (1*13) 

This  is  shown  by  curve  ■  c  in  figure  1 

Constant  strain  rote  is  specified  by  ^ 

■  c  t  H(t)  and  -  c  \ 


From  eq.  (1.1 2  ) 
^  (t)  -  L"1 


c*  1 

tllm  c~2 
P 


3  K  (|  +  p)  -  2  Gp 
2[  3  K(  i  +  p)  +  Gp)  p 


(1.13)' 
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1  0  1  r  1  3  K  t  i 

?‘2Ttl  1  *  exp{m  JTTZ  x  J 


(1.1b) 


I’his  is  shown  by  curve  ^  *  c  t  in  figure  1 


But  <rn 

I'hen  eq.  (1.13)  becones 


Jk  us  _  1  .  /  3  K  -  1  G  ,,  __  ,  3  K  t 

s  (t)  ? +  (  rmr  ■ 1]  exp  (_  rrrc- 


T/n 


and  eq.  (l.lh)  becones 


>)  (0  •  \  -  ^ 


/  juduj 

n nr 


l/n 


j  1 2 3 * 5 6 * 8  -  exp  (  - 


§T 


which  gives  the  ratio  cf  lateral  tc  longitudinal  strain  in  uni-axial 
stressing  in  t ems  cf  viscosity  and  fluid  velocity. 

The  paper  has  proposed  a  generalization  of  the  Navier-Stokes 
equation  for  non-Newtonian  fluids  in  which  the  stress  is  propor¬ 
tional  to  a  power  n  cf  the  rate  cf  strain,  i’he  ratio  of  the  trans¬ 
verse  to  the  longitudinal  strain  is  not  a  constant  but  a  function 
of  tia*;  or  cf  strain  and  depends  strongly  upon  the  visco-elastic 
response  of  the  medium  as  well  as  on  the  loading  conditions. 
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